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Abstract 

We propose a parametrization for two-body nonleptonic B meson decays, in which the various topologies 
of amplitudes are counted in terms of powers of the Wolfenstein parameter A ~ 0.22. The weak phases 
and the amplitudes are determined by comparing this parametrization with available measurements. It is 
possible to obtain the phase 03 from the B — > Ktt data up to theoretical uncertainty of O(A^) ~ 5%. The 
recently measured B^ — » 7r°7r° branching ratio implies a large color-suppressed or penguin amplitude, and 
that the extraction of the phase (j)2 from the B tttt data may suffer theoretical uncertainty more than the 
expected one, O(A^) ~ 5%. 



* Email: charngOphys. sinica. edu.tw 
^ Email: hnli@phys. sinica. edu.tw 



One of the major missions in B physics is to determine the weak phases in the Kobayashi-Maskawa ansatz 
for CP violation . The phase (j>i can be extracted from the CP asymmetry in the B J/ipKs decays in an 
almost model-independent way, which arises from the B-B mixing. The application of the isospin symmetry 
to the S — > TTTT decays j2j and to the B ^ pir decays |3j has been considered as giving a model-independent 
determination of the phase 02- However, this strategy in fact suffers the theoretical uncertainty from the 
electroweak penguin, which is expected to be about 5-10% . The phase <f>^ can be extracted in a theoretically 
clean way from the modes involving only tree amplitudes, such as i? ^ ttD 4 and B KD The difficulty 

is that one of the modes, such as i?° — > tt^ or B+ K^D^ , has a very small branching ratio and is not 
experimentally feasible [?]• The alternative modes B K*D |8| and B^ DgD W improve the feasibility 
only a bit. It has been pointed out that the B^ —> K^{D^ /) and B^ K^{D^ — > /) amplitudes, with 
l)" / being a doubly-Cabibbo suppressed decay, exhibits a strong interference [7| IIUL [TT] . For this strategy, 
the strong phase difference between f and — > / is a necessary input. Another possibility is to measure 

the B D*V decays for the vector meson V = p, K, • • •, since an angular analysis involves many observables, 
which are sufficient for extracting 03 model-independently |12) . 

Instead of resorting to theoretically clean modes, which are usually experimentally difficult, one considers 
the modes with higher feasibility and tries to constrain the decay amplitudes and the weak phases. The 
problem is that available measurements are usually insufficient to make the constraint, and theoretical inputs are 
unavoidable. For example, one adopts the (imaginary) tree-over-penguin ratio obtained from the perturbative 
QCD (PQCD) formalism 03 d [TSl UHl [I7| or from the QCD-improved factorization (QCDF) (TH|, so that 
the phase 02 can be extracted from the CP asymmetries of the i?^ — > tt+tt" decays. One may also employ 
symmetries to relate the amplitudes of the relevant modes, such as SU (3) ^Hl Etnd t/-spin (201) in order to reduce 
the number of free parameters. However, the theoretical calculations are subject to subleading corrections, and 
the symmetry relations are broken with unknown symmetry breaking effects. For these strategies to work, the 
theoretical uncertainty must be under control. 

In this paper we shall propose counting rules for the various topologies of amplitudes in two-body non- 
leptonic B meson decays in terms of powers of the Wolfenstein parameter A ^ 0.22 [32j. The relative importance 
among the topological amplitudes has been known from some physical principles: helicity suppression (color 
transparency) implies that tree annihilation (nonfactorizable) contributions are smaller than leading factoriz- 
able emission contributions. Here we shall assign an explicit power of A to each topology, such that the relative 
importance becomes quantitative. This assignment is supported by the known QCD theories |16[ 1181 1^ I24| . 
and differs from that assumed in 22 . We drop the topologies with higher powers of A until the number of free 
parameters are equal to the number of available measurements. The weak phases and the decay amplitudes can 
then be solved by comparing the resultant parametrization with experimental data. Afterwards, it should be 
examined whether the solved amplitudes obey the power counting rules. If they do, the extracted weak phases 
suffer only the theoretical uncertainty from the neglected topologies. If not, the inconsistency could be regarded 
as a warning to QCD theories for two-body nonleptonic B meson decays. For example, the long-distance rescat- 
tering effect has been neglected in PQCD and in QCDF. If this effect is important, the hierarchy among the 
various topological amplitudes will be destroyed |2SI ■ The comparison of our parametrization with data can tell 
whether the above assumption is reliable |2(ij . 

As shown below, dropping the electroweak penguin amplitude, the phase 02 can be extracted from the 
B ~+ TTTT data. In principle, the theoretical uncertainly of the ignored amplitudes is around O(A^) ~ 5%, the 
same as in the extraction based on the isospin symmetry 2 . Similarly, the phase 03 can be best determined 
from the B Ktt data up to the uncertainty from the neglect of the 0{X^) ^ 1% tree annihilation and 
color-suppressed electroweak amplitudes. Note that the determination of the phase 0i from the B J/^pK^*^ 
decays also bears about 1% theoretical uncertainty. Certainly, a CP asymmetry is an 0(A) quantity itself. 
Precisely speaking, the above determination of 02 and 03, involving the data of CP asymmetries, in fact carries 
the uncertainly of 0(A) ~ 20% and O(A^) ~ 5%, respectively. Because the B tttt, Ktt measurements are 
not yet complete, we shall drop more topologies in order to match the currently available data. In this simple 
demonstration, we observe that the amplitudes solved from the B Kn data more or less obey the hierarchy 
in A. That is, an almost model-independent determination of 03 is promising. The solution from the B tttt 
analysis is, unfortunately, not consistent with the power counting rules, indicating that the extraction of 02 may 
suffer theoretical uncertainty larger than stated above. Hence, our work casts a doubt to the strategy based on 
the isospin symmetry |5] and gives a warning to the QCD calculations of the B ^ tttt modes |17l 1181 1?7) . 
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We start with the B Kn decays. The branching ratio of a two-body nonleptonic B meson decay is written 



B{B ^ M1M2) = 



TB 



-\A{B^MiM2) 



(1) 



where the hght-meson masses mvr and niK have been neglected, and the B meson mass and the B meson lifetimes 
take the values ms = 5.28 GeV, tb± = 1.674 x 10^^^ s, Tgo = 1.542 x 10^^^ s. The effective Hamiltonian for 
the flavor-changing b ^ s transition is [2HI , 
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Ci{fi)0[''\fi) + C2{fi)0i''\fi)+J2c.{fi)0,{fi) , (2) 
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with the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements V and the operators, 



o['^^ = {siqj)v-A{qjbi)v-A , 
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q g 



Or = -isibi)v-A^eq{qjqj)v+A , Os ^ -{sibj)v~A'^eq{qjqi)v+A 
q q 

3 



O9 = :^{sibi)v-A^eq{qjqj}v-A , Oio = -{sibj)v^A^eq{qjqt) 



V-A 



(3) 



i, j being the color indices. For the characteristic scale /i ^ rribA ~ 1.5 GeV involved in two-body B meson 
decays A = 7715 — 77ib being the B meson and b quark mass difference, the values of the Wilson coefficients 
are 
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The above characteristic scale has been confirmed by the dynamical penguin enhancement exhibited in the 
B VP data 125. The Wolfenstein parametrization for the CKM matrix is given by 
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(5) 



with the parameters A = 0.2196 ± 0.0023, A = 0.819 ± 0.035, and = ^fp'- + 77^ 



= . /.2 ^ „2 _ 0.41 ± 0.07 |2i. Note that 
the product AR^ ~ 0.3 should be regarded as being of 0(A), and that jKbl is in fact 0(A*). The phases 01 and 
(/)3 are defined via Vtd = |Vtd| exp(— i(/)i) and Kb = |Kb| exp(— 1^3), respectively. 



Considering all possible topologies of amplitudes, the B — > Kn decay amplitudes are given by 
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which satisfy the quadrangle relation, 



A{B+ K\+) + V2A{B+ K+t:") = ^(B^ ^ A'+tt") + \/2A(B3 ^ K\°) . (10) 

The amplitude P^w [Penn Pew) comes from the color- allowed (color-suppressed, annihilation) topology through 
the electroweak penguin operators. The amplitude P includes the emission and annihilation topologies through 
both the QCD and electroweak penguins: 

P = Pqcd + euP:^ + euP^^ , (11) 

with the u quark charge e„ = 2/3. The amplitude T (C, T°) comes from the color-allowed (color-suppressed, 
annihilation) topology through the tree operators. The penguin contributions from the c quark loop can be 
included using the relation V*gVcb = —V*gVub — Ks^tfci ^-^id the expressions in Eqs. ((T7|) - (pn|l remain unchanged. 

It has been shown in PQCD that a nonfactorizable amplitude M"-^, a factorizable annihilation amplitude 
P{v-A) fro™ the {y — A){y ~ A) current, and a factorizable annihilation amplitude F^y_^j^-^ from the {V — 
A){y H- A) current are suppressed, compared to the leading factorizable emission amplitude F'' ^ by the factors 

ofism, 



, ruB 
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^ ^{v~A) ^ Aqcd _ ^2 ^{V+A) ^ 2mo ^ ^0 

' F'^ TTlB ' F'^ TTlB ' 



respectively, where ttiq is the chiral enhancement scale, and the CKM matrix elements and the Wilson coefficients 
are excluded. We list the power counting rules for the Wilson coefRcients in Eq. 

0(1) :ai , 

0(A) : a2 , 1/iVc , 

O(A^) : C4 , Cg , 04 , flg , 

C3 , C5 , Cg , 03 , 05 , ag , 
Cio , 

C7 , Cs , 07 , as , aio , (13) 



with ai = C2 + Ci/N„ a2^Ci + C2/N,, a, ^ G, + C,+i/N, for z = 3, 5, 7, 9, and a, = + Q-i/N, for 
i = 4, 6, 8, 10. 

According to Eqs. l(T^ and l(T^. we assign the powers of A to the following ratios of the various topological 
amplitudes: 
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(14) 



For the latter three ratios, we present the power counting rules derived from both the factorizable and nonfac- 
torizable contributions, which are of the same order of magnitude. Compared to the power counting rules in 
22] based on the conventional scale /i ^ mf,, P^^/P is down by one more power of A due to oio ^ O(A^) in 
PQCD. 

Whether a factorizable amplitude or a nonfactorizable amplitude is important depends on the decay modes. 
In the B Ktt case, C mainly comes from the factorizable color-suppressed diagrams, since there is a strong 
cancellation between a pair of nonfactorizable diagrams. The factorizable and nonfactorizable annihilation 
contributions to T", P^^, and P"^ are of the same order of magnitude as shown in Eq. lfn|l . In the B Dw 
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decays, C, being of the same order of magnitude as T, mainly comes from the nonfactorizable color-suppressed 
diagrams, since the above cancellation does not exist |15[ \'61\ . For T" in the B — > Dtt case, the nonfactorizable 
diagrams dominate, because of 



Employing the reparametrizations, 
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we arrive at the most general parametrization of the B Kt: decay amplitudes. 
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There are totally 6 independent amplitudes, namely, 11 unknowns, because an overall phase can always be 
removed. Hence, we choose the amplitude P as a positive real value. Plus the weak phase (/)3, the 12 unknowns 
are definitely more than the 9 experimental inputs: the branching ratios and the direct CP asymmetries of 
the four modes, and the mixing-induced CP asymmetry of the B^ — s- X^tt'^ mode. Dropping the 0{\^) terms, 
T"" /P and P^^/P, we have 8 unknowns. Then the data of the direct CP asymmetry in the P+ — > K^tt^ decays 
should be excluded for consistency. Hence, we have 8 experimental inputs, and thus all unknowns can be solved 
exactly assuming the phase 0i is already known from the measurement of the mixing-induced CP asymmetry 
in the B J/ipK^*"^ modes. The determination of 03 from this parametrization is then accurate up to the 
theoretical uncertainty of O(A^) ^ 5%. 

We emphasize the consequence from the different power counting rules in |22| and in this work: the smaller 
Pg^ is crucial for claiming that the determination of (/)3 from the B — > Kit data is accurate up to 5% theoretical 
uncertainty. Following the counting rules in ^2^, both C and P^^ will be included at O(A^), such that the 10 
unknowns are more than the 9 available measurements. In this case we can not solve for C and P^^ exactly, and 
have to rely on symmetry relations to reduce the number of unknowns. It is then difficult to estimate the involved 
theoretical uncertainty. Using the counting rules in Eq. (|14ll . which are supported by the PQCD calculation, we 
include only C at O(A^), and the number of unknowns can be equal to the number of measurements. Solving 
for C, and assuring that the solution obeys our counting rule as a self-consistency check, the uncertainty from 
the neglected topologies is under control. 

The measurement of the time-dependent asymmetry in the B^ — > Kstt'^ decay still suffers a large error. 
To demonstrate our method, we reduce the number of unknowns by further dropping the O(A^) terms, C/P, 
arriving at 
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(21) 
(22) 

(23) 

(24) 



where St and Sew denote the strong phases of T and Pew, respectively. The B — + Ktt decay amplitudes in 
Eqs. (|21|) - (|24|l are the expansion up to the power of A, at which the determination of 03 suffers the theoretical 
uncertainty of 0(A) ~ 20%. 
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Figure 1: The allowed range of T/P determined from the B Ktt data. 

Re [Pew/P] 




Figure 2: The allowed range of Pew/P determined from the B — > Ktt data. 



We shall solve for the 6 unknowns: 
experimental inputs , 



\T\, 03, and the strong phases 5ew and 5t, from the 6 
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A{B^ - 
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* K°Tr^) = (11.5 ± 1.7) X IQ-s 
i^±7r±) = -(10.2 ±5.0)% , 

> K^tt") = -(9.0 ±9.0)% . 



(25) 



and B^ 



K^TT^ modes indeed have very small direct CP asymmetries, consistent with 



The bounds on the various amplitudes and phases can be derived 



The B± iC°7r± 

the parametrization in Eqs. H21|l - H24|) . 
unambiguously from Eq. (|25|l . 

The allowed ranges of the ratios T/ P and P^w/P are exhibited in Figs. Handle respectively. The prescription 
for deriving the two figures is briefly explained below. The data for each branching ratio and for each CP 
asymmetry are expressed as a set, whose elements are the central value with +lx error bar, Ox error bar, and 
— Ix error bar. For a combination of the element from each set, we solve the coupled equations, and the solution 
is represented by a dot in the figure. Scanning all the combinations, we obtain the ranges in the figures. The 
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central values of the solutions are 
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0.23 , 5t = - 

= 0.50, Sen, 



13° 



(26) 



The above resuh of T/P is in agreement with the PQCD prediction, T/P - 0.20 exp(-27°z) [HEll El ES] , 
while the central values of \Pew\/P and of Sew differ from the PQCD prediction, |Pe^|/P ^ 0.2 and Sew ~ St, 
respectively. The latter PQCD prediction is consistent with the almost model-independent relation between the 
electroweak penguin and tree amplitudes obtained in |25[I36| . The ratio \Pew\/P — 0.5 and the nearly 90° phase 
between Pew and P in the above fit have been speculated in I37L EHl • We also derive the allowed ranges 
0.06 < \T\/P < 0.72 and 0.22 < \Pew\/P < 0.70, implying that the extracted ratios \T\/P and \Pew\/P deviate 
a bit from the power counting rules in Eq. H14|l. Hence, the B —>■ Kn data are indeed puzzling, especially from 
the viewpoint of the dramatically different strong phases Sew and St shown in Figs. ^ and 13 Because of the 
large central values of \Pew\/P and of Sew, a strong hint of new physics has been claimed in jSTI 1381 139| . A 
more convincing examination of the self-consistency can be made by solving for the amplitude C, when more 
complete data are available. At last, the central value and the allowed range of the phase 03 are given by 



102° 



26° < 03 < 151° , 



(27) 



respectively, with the theoretical uncertainty of about 20%. 

We emphasize that our fitting differs from the global fitting based on the QCDF approach I41|. For 
example, the penguin contributions have been split into the factorizable type depending on a transition form 
factor, the nonfactorizable type depending on the imaginary infrared cutoff pn for an end-point singularity, 
and the annihilation type depending on the imaginary infrared cutoff pA in QCDF. Taking into account only 
the B PP modes, such as i? ^ Ktt and tttt, the fitting result of the phase 03 ^ 110° is close to that 
extracted in this work. Our method also differs from those based on the isospin relations 02]) with which some 
combinations of the B Kn branching ratios can be described by the functions of the parameters Pew/T, T/P 
and the relative strong phases. The SU{3) flavor symmetry is then employed to fix Pew/T and T/P. Finally, 
only the strong phases and the weak phase 03 are treated as unknowns, and determined by the data. The 
conclusion is similar: the B Kt: data favor 03 > 90°. Our approach does not rely on the SU{3) symmetry, 
and the ratios Pew/T and T/P are treated as unknowns. Including the B — > VP modes in the QCDF fitting, 
the value of 03 could be smaller than 90° @T] . Using the parametrization for the the B VP modes based on 
SU(3) flavor symmetry, an phase 03 < 90° was also obtained ^]. In a forthcoming paper we shall apply our 
parametrization to the B — > VP modes, and make a comparison with the above works. 

Next we apply our method to the B ^ tttt decays. Considering all possible topologies of amplitudes, their 
decay amplitudes are given by 
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which satisfy the triangle relation, 

V2A{B+ 7r+7r°) 
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In the above expressions the amplitude P has been defined in Eq. Hll|) . and the annihilation contribution 
comes only from the nonfactorizable diagrams. Based on Eqs. (|12|l . (|13|l and 1)14(1 . we assign the power counting 
rules to the following ratios of the topological amplitudes: 
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Employing the reparametrizations, 
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the most general parametrizations of the B —f tttt deeay amplitudes are written as 
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There are 4 independent amplitudes, namely, 8 parameters including the phase 02, which are more than the 
available measurements. Neglecting the O(A^) term, P^w/T^ the resultant expressions are the same as in |22| : 



^/2A{B+ - 
V2^(i?° 



7r+^") = 



ttV") 



1^1 .rSc 



-T\l + ye 



\P\ 
T 



-T ( 1 + H-Le^'^^e*''^ 



\P\ 
T 



\C\ 



,iSc 



(37) 
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for which we have 6 unknowns T, |C|, |P|, 6c, Sp and 02- Similarly, we have removed the strong phase of T, 
and assumed it to be real and positive. 

In this case we have to exclude the data of the direct CP asymmetry in the i?+ ^ 7r+7r° decay, and 6 
experimental inputs are relevant: the three CP-averaged branching ratios, the direct and mixing-induced CP 
asymmetries in B^ tt+tt", and the direct CP asymmetry in B'^ 7r°7r°. At this level of accuracy, our 
parametrization is equivalent to that based on the isospin triangle [2 I47| , in which the electroweak penguin 
contribution to the B^ — > tt^tt'^ decay is also ignored. We mention that the electroweak penguin amplitude has 
been included in the isospin analysis of the B — > tttt decays, and that the CP asymmetry in the B^ — > tt^tt^ 
modes still vanishes After extracting (j)2 from the i? — > tttt data and (j)^ from the B — > Kn data, we can 
check whether they, together with from the B J/^pK^*^ data, satisfy the unitarity constraint, when the 
data precision improves. 

The time-dependent CP asymmetry of the B^ tt+tt^ mode is expressed as 
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where the direct asymmetry C^rn and the mixing-induced asymmetry S'jr^ are defined by 
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Atttt — e 



2id 



1 + e-"t'^P/T 
1 + e^'i'^P/T 



(40) 



(41) 



(42) 



7 



Re [P/T] 



Figure 3: The allowed range of P/T determined from the B ^ tttt data. 
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Figure 4: The allowed range of C/T determined from the B ^ tttt data. 
The data are summarized as j46| . 

>7r±7r°) = (5.2 ±0.8) x 10"^ , 
7r±7r=F) = (4.6 ±0.4) x 10"^ , 
ttOtt^) = (1.97 ±0.47) x 10"^ , 

(38 ± 16)% , = -(58 ± 20)% . (43) 

Since the data of the direct CP asymmetry in the B'^ tt^tt" mode is not yet available, we shall assign a 
plausible range to it, 

A{B"^ ^ vr^vr") = (-50 - +50)% . (44) 

The central values of the measured B^ tt^tt^ and B^ tt^tt^ branching ratios are close to each other, 
implying that either C is large and constructive in order to enhance the B^ n^n^ modes, or P is large and 
destructive (after including the weak phase 02 ) in order to suppress the 5° tt^tt^ modes [^T]. In either case 
the 5° tt'^tt'^ branching ratio exceeds the expected order of magnitude, 0(10"''). There exist four solutions 
associated with each set of data input: two solutions correspond to the large C and P cases, and the other two 
are the reflections of the first two with respect to the B^ tt^tt'^ side of the isospin triangle. Note that the 
relations of the phase (/)2 to the measured quantities have been given in 03 without numerical results. Here we 
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shall not present the central values of the solutions, because the central values of the experimental data of the 
7r°7r" direct CP asymmetry are not yet available. 
The ranges of P/T and C/T, shown in Figs.|niand01 respectively, collect all allowed solutions. These ranges 
indicate that the hierarchy in Eq. H32|) is not satisfied, since both \P\/T and \C\/T can be as large as 1, much 
greater than 0(A) ^ 0.22. There is then no reason for believing that the effect of the electroweak penguin would 
be as small as O(A^) ~ 5% according to the relation between P^w and T |25[ Our analysis implies that the 
extraction of 02 from the B —^ tttt data based on the isospin symmetry may suffer the theoretical uncertainty 
more than expected. It also casts a doubt to the PQCD (also QCDF) calculation of the B ^ tttt decays. To 
complete our numerical study, we present the allowed range of (j>2 corresponding to the data in Eq. (|43|l . 

51° < 02 < 176° . (45) 

As explained above, the theoretical uncertainty associated with the above range may not be under control. 

When the data become more precise, and when the data of more CP asymmetries, such as the mixing- 
induced CP asymmetry in the B^ Kstt'^ mode, are available, the allowed range will shrink, and the theoretical 
uncertainty can reduce. Our method then tells whether the B — > Ktt data indicate a solid signal of new physics. 
Besides, our parametrization extends straightforwardly to the other relevant modes, such as B _fC*7r, from 
which the phase 03 can also be extracted j4Q^. Considering the overlap of the extractions from different modes, 
the allowed ranges of the decay amplitudes and of 03 will shrink too. An evaluation of the next-to-leading-order 
corrections to the B tttt decays in the PQCD framework is now in progress, whose result will clarify whether 
the large |P| or |C| is understandable. If not, new dynamics, such as the rescattering effect, might be important. 
The B — !■ TTTT decays and the extraction of the phase 02 then demand more theoretical effort. 



We thank R. Fleischer, A. Hocker, Y.Y. Keum, H. Lacker and D. London for useful discussions. This work 
was supported by the National Science Council of R.O.C. under the Grant No. NSC-92-2112-M-001-030 and 
by the National Center for Theoretical Sciences of R.O.C 

References 

[1] M. Kobayashi and T. Maskawa, Prog. Th. Phys. 49, 652 (1973). 
[2] M. Gronau and D. London, Phys. Rev. Lett. 65, 3381 (1990). 

[3] E. Snyder and H.R. Quinn, Phys. Rev. D 48, 2139 (1993); H.R. Quinn and J.P. Silva, Phys. Rev. D 62, 
054002 (2000). 

[4] R. Aleksan, I. Dunietz, B. Kayser, and F. Le Diberder, Nucl. Phys. B361, 141 (1991); I. Dunietz, Phys. 
Lett. B 427, 179 (1998) 

[5] M. Gronau and D. Wyler, Phys. Lett B 265, 172 (1991); M. Gronau and D. London, Phys. Lett B 
253, 483 (1991); A.I. Sanda, hep- ph/0108031| M. Gronau, Phys. Lett. B 557, 198 (2003); R. Fleischer, 
'hep-ph/0304027, 

[6] I. Dunietz, Z. Phys. C 56, 129 (1992). 

[7] D. Atwood, I. Dunietz and A. Soni, Phys. Rev. Lett. 78, 3257 (1997); Phys. Rev. D 63, 036005 (2001). 

[8] I. Dunietz, Phys. Lett. B 270, 75 (1991). 

[9] R. Fleischer and D. Wyler, Phys. Rev. D 62, 057503 (2000). 



[10] B. Kayser and D. London, (hep-ph/9905561 D. Atwood and A. Soni, 1hep-ph/0206045| |hep-ph/0212071 
|hep-ph/0304085, 

[11] Y. Grossman, Z. Ligeti, and A. Soffer, Phys. Rev. D 67, 071301 (2003); R. Fleischer, Nucl. Phys. B659, 
321 (2003). 



9 



[12] D. London, N. Sinha and R. Sinha, Phys. Rev. Lett. 85, 1807 (2000); N. Sinha and R. Sinha, Phys. Rev. 
Lett. 80, 3706 (1998). 

[13] H-n. Li and H.L. Yu, Phys. Rev. Lett. 74, 4388 (1995); Phys. Lett B 353, 301 (1995); Phys. Rev. D 53, 
2480 (1996). 

[14] C.H. Chang and H-n. Li, Phys. Rev. D 55, 5577 (1997). 
[15] T.W. Yeh and H-n. Li, Phys. Rev. D 56, 1615 (1997). 

[16] Y.Y. Keum, H-n. Li and A.L Sanda, Phys Lett. B 504, 6 (2001); Phys. Rev. D 63, 054008 (2001). 
[17] CD. Lii, K. Ukai, and M.Z. Yang, Phys. Rev. D 63, 074009 (2001). 

[18] M. Beneke, G. BuchaUa, M. Neubert and C.T. Sachrajda, Phys. Rev. Lett. 83, 1914 (1999); NucL Phys. 
B591, 313 (2000); NucL Phys. B606, 245 (2001). 

[19] M. Gronau, J.L. Rosner, and D. London, Phys. Rev. Lett. 73, 21 (1994). 

[20] R. Fleischer, Phys. Lett. B 459, 306 (1999); M. Gronau and J.L. Rosner, Phys. Lett. B 482, 71 (2000). 
[21] L.L. Chau, H.Y. Cheng, and B. Tseng, Phys. Rev. D 43, 2176 (1991). 

[22] M. Gronau, O.F. Hernandez, D. London, and J.L. Rosner, Phys. Rev. D 50 4529, (1994); Phys. Rev. D 52 
6356, (1995); Phys. Rev. D 52 6374, (1995). 

[23] Y.Y. Keum and H-n. Li, Phys. Rev. D63, 074006 (2001). 

[24] C.H. Chen, Y.Y. Keum, and H-n. Li, Phys. Rev. D 64, 112002 (2001); Phys. Rev. D 66, 054013 (2002). 

[25] M. Neubert, Phys. Lett. B 424, 152 (1998); J.M. Gerard and J. Weyers, Euo. Phys. J. C 7, 1 (1999); C. 
Smith, lhep-ph/0309062, 

[26] C.H. Chen and H-n. Li, Phys. Rev. D 63, 014003 (2001). 

[27] M. Beneke and M. Neubert, hep-ph/0308039", 

[28] G. BuchaUa, A. J. Buras and M. E. Lautenbacher, Review of Modern Physics, 68, 1125 (1996). 

[29] Review of Particle Physics, Eur. Phys. J. C 3, 1 (1998). 

[30] H-n. Li and K. Ukai, Phys. Lett. B 555 197, (2003). 



[31] Y.Y. Keum, T. Kurimoto, H-n. Li, CD. Lu, and A.I. Sanda, |hep^p h/0305335' 
[32] Y.Y. Keum and A.I. Sanda, |hep-ph/0306004| 



[33] BaBar Collaboration, B. Aubert et al., hep-ex/0207055 



[34] K. Ukai and A.I. Sanda, Prog. Theor. Phys. 107, 421 (2002). 



[35] Y.Y. Keum and A.I. Sanda, Phys. Rev. D 67, 054009 (2003); Y.Y. Keum, |heM)h/0209208 
[36] M. Neubert and J.L. Rosner, Phys. Rev. Lett. 81, 5076 (1998). 
[37] M. Gronau and J.L. Rosner, |hep-ph/030709'5l 



[38] A.J. Buras, R. Fleischer, S. Recksiegel, and F. Schwab, |hep-ph/0309012 
[39] T. Yoshikawa, Phys. Rev. D 68, 054023 (2003). 



[40] W.M. Sun, hep-ph/0307212 



10 



[41] D. Du, J. Sun, D. Yang, and G. Zhu, Phys. Rev. D 67, 014023 (2003); R. Aleksan, P.P. Giraud, V. Mornas, 
O. Pne, and A.S. Safir, Phys. Rev. D 67, 094019 (2003). 

[42] M. Beneke, hep-ph/0207228 

[43] R. Fleischer, | hep-ph/0306270| and references therein. 
[44] C.W. Chiang et d., |hep-ph/0307395 ) 

[45] M. Gronau, D. Pirjol, and T.M. Yan, Phys. Rev. D 60, 034021 (1999). 

[46] H. Jawahery, talk presented at the XXI International Symposium on Lepton and Photon Interactions at 
High Energies, Fermilab, USA, Aug. 2003; BaBar Collaboration, B. Aubert et ai, h ep-ex/0208012 , 

[47] H.J. Lipkin, Y. Nir, H.R. Quinn, and A.E. Snyder, Phys. Rev. D 44, 1454 (1991). 



Z.Z. Xing, hep-ph/0308225 



11 



